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Abstract: In a large class of extra-dimensional models, a scalar degree of freedom known as 
the radion is long-lived, or even stable, on cosmological scales. In this paper we investigate 
the impact of radionactivity on the evolution of the universe. We demonstrate that whether 
the radion overdoses the universe, constitutes the dark matter, is the inflaton, the curvaton, 
or does not play any role in cosmology, depends crucially on the ratio between the energy 
densities stored in the radion and in the inflaton at the time of inflation. We discuss the 
general difficulties reconciling models with low compactification scale (i.e., TeV scale) with 
the simple picture of inflation. 
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1. Introduction 



Scalar fields are believed to play important roles in both particle physics and cosmology. The 
crucial aspect of scalar fields is that they determine the properties of the vacuum without 
relinquishing the key property: invariance under spacetime symmetries. In spite of these very 
desirable properties, scalar fields fit rather uneasily within quantum field theories. Scalar field 
masses are extremely sensitive to quantum corrections. They exhibit quadratic divergences, 
which suggests that their natural mass scale should be very large, and implies that their 
detailed properties should be controlled by ultraviolet (UV) physics beyond our ability to 
infer from current experiments. The large masses natural for scalar fields tend to spoil their 
potential utility in both particle physics and cosmology. Idealy, we should be led to consider 
theories with some mechanism to protect the scalar mass from large corrections. 

In this article we explore one generic feature of theories with extra spatial dimensions. 
Such theories contain a scalar field, the radion, which plays the role of the geometrical mod- 
ulus that determines the size of the compactified dimensions. From the higher dimensional 
perspective, the radion is part of the enlarged graviton tensor, and at distances smaller than 
the size of the extra dimension it is protected from receiving a large mass by the same general 
covariance which protects the graviton mass. Thus, unlike generic scalar fields, it can be 
naturally light. Light scalar degrees of freedom are often invoked in cosmology: inflatons, 
curvatons, modulons (fields controlling coupling strengths), dark matter, etc. It is interesting 
to explore the possibility that the radion may be one of the much-needed light scalar degrees 
of freedom, and to see if it can play an interesting role in cosmology. 

Extra-dimensional models with a TeV compactification scale [|J have received revived 
interest lately. There are many ways in which extra dimensions can attempt to explain the 
outstanding mysteries in particle physics [§, [3|, |Q. For example, there is the prospect of 
dark matter consisting of the first Kaluza-Klein (KK) mode of the hypercharge gauge boson 
(B 1 ) H in a sub-class of these models referred to as universal extra dimensions (UED) ||. 
However, the cosmology of models with TeV size extra dimensions remains largely unexplored. 
Our concern in this paper is to investigate, for all possible compactification scales, constraints 
from radion cosmology on these models, and to a larger extent, on any models with flat extra 
dimensions. Different recent attempts discussing the role of the radion during inflation can 
be found in |^] . 

The destiny of a radionactive universe depends on two parameters, the compactification 
scale M c (naturally L , the physical size of the extra dimensions) and the ratio r = M c /Mj 
where Mj = Vj^ is the scale of the inflaton potential. One question to ask concerns the epoch 
at which the extra dimensions were stabilized compared to the epoch of inflation. When did 
compactification take place relative to inflation? A first possibility is that inflation took place 
before compactification. This would require a post-inflationary exponential contraction of 
extra dimensions, and seems counterintuitive. Further, it would require that the inflation 
scale M/ is much larger than the compactification scale M c . As is well known, the non- 
renormalizability of extra-dimensional theories implies that they have a low cut-off scale A, 
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beyond which the theory is strongly coupled and loses predictivity. 1 

The second possibility stabilizes the extra dimensions, and then inflates the four large 
ones much later. In this case, geometrical moduli such as the radion are active during in- 
flation. These are scalar fields with Planck-suppressed couplings, and may contribute to the 
inflaton potential. In fact, they may be responsible for triggering inflation, for all possible 
compactification scales || . Even if they do not have the properties required to trigger infla- 
tion, they may still play a crucial role in the generation of cosmological perturbations and/or 
make up the dark matter. As we shall see, these scalar fields are typically light, with masses 
protected by higher dimensional general covariance. It is well known that when the mass 
is smaller than the Hubble rate during inflation, the field is frozen during inflation. After 
inflation, the energy density stored in these fields can easily overdose the universe. This is 
the essence of the moduli problem. Another well known fact is that during inflation, scalar 
fields receive corrections to the their mass of order H from non renormalizable operators of 
the type Vi4> 2 /m 2 pi (where <fi is the modulus in question). Scalars with masses larger than ap- 
proximately H are damped out to the minimum of their potential during inflation. However, 
their energy is not necessarily inflated away because the position of the minimum is shifted 
by the same type of operators which are responsible for shifting the masses. As a result, 
energy stored in the radion fields typically survives inflation and can have dramatic effects on 
post-inflationary cosmology. We can actually predict m 2 /H 2 during inflation as a function 
of r. We show that if r < 1, models with a compactification scale M c < 10 9 GeV can easily 
be ruled out (at least within a given inflationary paradigm) because of the overclosure of the 
universe by radion oscillations. This means that inflation must take place at a scale as low 
as or lower than the compactification scale. 

Inflation taking place after compactification appears problematic in models with a low 
compactification scale, i.e., M c ~ TeV. This implies that the highest energy scale of the 
universe when inflation starts is Mj < M c . From the COBE constraint 5h ~ 1CP 5 [ 10 1 and 



by identifying the energy density fluctuations with primordial quantum fluctuations of the 
inflaton, one obtains: 



e 

,2 nr. i nr. \z 



^ ...16 



10 iD GeV, (1.1) 



where e = m pi {Vi' /Vi) /2 is the slow-roll parameter. This data disfavors models with low- 
scale inflation. Equation ( |1.1| ) tells us that if the primordial density perturbations originate 
from the inflaton, inflation at the TeV scale requires a slow-roll parameter of the order of 
10~ 48 , a tiny value which is difficult to motivate naturally. However, it has recently been hy- 
pothesized that structures in the universe may not necessarily follow from primordial quantum 
fluctuations in the inflaton, but could instead be due to another scalar field, referred to as the 
curvaton | |1 lfl . The curvaton scenario has the advantage that the constraints on the inflaton 
potential are relaxed, though at the cost of introducing another scalar field whose mass must 

1 For theories in which gauge fields propagate in the bulk, naive analysis suggests A ~ 100M C (see however 
Q which suggests A ~ M c ) if there are 5 dimensions, and A is closer to M c for theories with more dimensions. 
The scale at which gravitational interactions become strong is typically about M| ~ M c m 2 Pl 3> A 3 in 5d. 
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somehow be protected. Furthermore, the constraints are not relaxed to the point that TeV 
inflation becomes natural. Indeed, the COBE constraint in the curvaton scenario translates 
into 

« 10" 4 , (1.2) 

where and 0* are the Hubble constant and the amplitude of the curvaton field respectively 
at the time of horizon exit. Here fu is the fraction of the total energy density contained in the 
curvaton when it decays and should be of order unity for this mechanism to work. To satisfy 
Eq. (|l.2| ), an inflation scale at a TeV would require the small value 0* ~ 10~ 9 GeV <C Mj 
for the curvaton amplitude, which is hard to reconcile with the possibility that the curvaton 
dominates the total energy density at the time of its decay. An even more recent proposal 

suggests that cosmological perturbations could be due to primordial fluctuations in the 
field (which we will refer to as a modulon) whose VEV gives the coupling between the inflaton 
and the Standard Model fields and plays a crucial role in the reheating process. In this work 
we will investigate whether the radion could ever play the role of a curvaton or a modulon. 

A third alternative has compactification and inflation taking place at comparable scales. 
We assume that the relevant history of our universe starts at a time when all dimensions 
have comparable sizes. Inflation starts in all 4 + n dimensions (with a typical scale around 
or a few orders of magnitude below mpi). Later on, the extra dimensions stop inflating, for 
example when M c is at the TeV scale, and the rest of the 4 dimensions continue to inflate. 
One assumes there are some extra fields living in the extra dimensions responsible for the 
stabilization dynamics and allowing the extra dimensions to undergo only 15 — 20 e- folds of 
inflation while allowing the other dimensions to continue to inflate. In this case, the scale 
in the inflaton potential is still much higher than the compactification scale, Mj 3> M c , 
and there would be no need for a huge fine-tuning in the epsilon parameter or in 0* to 
get the right COBE normalization. As we will see in section [f| there is a serious moduli 
problem associated with this scenario. Indeed, when the ratio between the compactification 
scale and the inflaton scale is much smaller than one, we are in the regime where there is no 
damping at all of the radion field during inflation, and it will overdose the universe. Avoiding 
the overclosure problem typically requires Mj <C M c , which on the other hand disfavours 
low compactification scales from the COBE constraint, unless there is a new mechanism 
to generate density perturbations. Obviously, there are interesting cosmological challenges 
associated with TeV compactification scales. In this work, in addition to examining the 
difficulties associated with low compactification scales, we will be able to constrain even 
larger compactification scales, up to 10 12 GeV. This leads us to study of more traditional 
compactifications models, i.e., associated with smaller extra dimensions. 

In summary, as soon as compactification took place before inflation, we are presented with 
additional scalar fields during inflation, whose masses are protected from large corrections, and 
thus they typically survive inflation. Note that as the number of extra dimensions increases, 
the number of geometric moduli also increases. We concentrate on the volume-modulus, i.e., 
we focus on the 5d case for illustrative purposes. In principle, one may also consider more 
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extra dimensions, and those moduli which describe the shape of the compact manifold. 

Section [2] reviews the definition and properties of the radion field and contains a few 
remarks about the possibility that the radion is itself the inflaton. Section 3 stresses the 
consequences of the generic form of the radion coupling to any form of energy density and 
how the radion can be naturally induced to be an inflaton. In section |] we evaluate the 
energy density stored in the radion at the end of inflation and discuss its evolution when it 
begins to oscillate. For small enough M c , the radion is stable on the time-scale of the age of 
the universe, and could either account for the dark matter, or overdose the universe. This is 
explored in section |5[ In section ^, after studying the radion isocurvature perturbations, we 
discuss whether the radion can be a viable curvaton or modulon candidate, and can account 
for the observed density perturbations. Section 7 studies the constraints associated to the late 
decay (after BBN) of radions and KK gravitons. In section 8 we summarize the challenges 
associated with TeV compactification scales and conclude. 

2. Radion Properties 

In this section, we briefly review the crucial properties of the radion: its mass, and its coupling 
to bulk fields. When reducing the curvature term of the (d = d + ra)-dimensional action, we 
obtain terms of the form 

_ _L J # x y[tfb'R@ = j d d x^gV)^)e^ -^2 nid) + \^9i3 F U F ^ V 

- ^rr^ (^)(^) - ^^'(teX^s) (2.1) 

+ -LgV(d^ k )g kl (d»g u ) - -^g ii (d fl g ij )9 kl (&*g kl ) + •• 

where the general metric is given by 

ds 2 = (e s <v + Agij^A^ADdx^dx" + g ij (2KA i VL dxPdy> + 2KA{dx v dy i ) + g ij dy i dy j . (2.2) 

All fields are functions of both extra-dimensional coordinates yi and four-dimensional coor- 
dinates x^. However, the zero- mode fields are independent of the yi, and thus we may drop 
this dependence when describing the low-energy (zero-mode) physics. In five dimensions, 
the KK tower of the radion is eaten by the massive KK gravitons leaving only the radion 
zero mode; in more extra dimensions, some of the radion KK modes remain as well, k is 
Newton's constant in d dimensions and is related to Newton's constant in d (flat) dimensions 
(k = nip 2 where mpi is the reduced 4-d Planck mass) by k = k/\/Vn where V n is the volume 
of the compact space. Here j" and g^ are the determinants of the metrics g^ v and gij. 
In the right-hand side, fi and v indices go up and down with the g^ v metric and its inverse. 
To obtain a curvature term in d dimensions which is canonically normalized, we make the 
following rescaling: 

5 (n) e ^ s = l ; i. e . ; e s = l/J~g~W . (2.3) 
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In this case, Eq. fl2.1|) simplifies to 

-4 I ' d i x\[g~^K^ = fd d xJgV) + l ff W^ 5ij ^> (2.4) 

K J ' J \_ K 4 



We now apply these results to the 5-dimensional case, S 1 /Z2 with bulk fields. The 5d 
action can be expressed as 

£5 = {M|i? 5 - ^Tr[F MN F MN ] - MlvS) - X -M%d M ^d M ^ - M$,%$)j , (2.5) 

where g, M5, and R§ refer to the 5d gravitational quantities, F MN is the field strength of a 
bulk gauge field with 5d coupling 55, i/j represents a minimally coupled real bulk scalar with 
potential and is the radion potential. For simplicity, we have neglected the possibility 
that tp mixes with the Ricci scalar; we will return to this later. 

A convenient parametrization for the 5d metric can be written in the form fl~3[| : 

ds 2 = (e'^g^ + e 2 ^A^A v ) dx^dx v + 2e 2/ ^A fl dx»dy + e 2/ ^dy 2 , (2.6) 



where = 2kA^ is the gravi-photon field, which is odd under the orbifold and thus has no 
zero mode, g^ is the tensor graviton field, and 4> is the (dimensionless) gravi-scalar (radion). 
This form is a convenient starting point because it decouples the kinetic terms of the radion 
and the transverse tensor modes and produces a radion kinetic term of the canonical form. 

To determine the effective action at scales less than of order 1/L, we replace the 5d fields 
with their zero-mode components, and integrate out the extra dimension. In doing so, it is 
convenient to introduce a reference length L. Note that the physical size of extra dimensions 
is Lphy S = e^/ 3 L. The effective action for the zero modes is, 



£4 = Vg 



, MlL ~ ~ L 7,, Mf,L ~ _ 



(2.7) 



where now all quantities refer to the 4d (zero-mode) fields. From here it is simple to rescale 
fields to canonical normalization and define four-dimensional scalar potentials, 

VM) = LM^VS), W) = LM$%$). (2.8) 

In order for the low energy physics to appear four dimensional, the radion potential must 
provide a stable VEV for (p. We can adjust the reference quantity L such that this happens 
for ((/>) = 0. Clearly, unless there is fine-tuning in the radion potential, this will occur for 
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L ~ M c 1 . For this choice, L is in fact the physical size of the extra dimension, and we can 
identify M|L = m Pl and g\jL = g\. The radion couplings are thus, 

- le^^F^F^ - e-^^K/XVO, (2.9) 

where we have rescaled the kinetic terms for the bulk fields to canonical normalization as in 
Eq. 



2.1 Radion Mass 

In the absence of contributions from V^h, which we consider in the next section, the radion 
mass may be determined by expanding the potential V$ appearing in Eq. ( |2,7D to order (f> 2 
about its minimum, 



-MlL-V;m 2 + ■■■ = -MlL-V'^)^ 2 + ... (2.10) 

where primes refer to functional derivatives with respect to <fi. Without fine tuning, one 
expects L ~ M~ , and all derivatives of to be order unity. Thus, the radion mass is 
naturally 

/-Ml , 
m r ~ V3 — -. (2.11) 
mpi 

Note that loop corrections will not destabilize the radion mass, as they are of order M 2 /mpi. 
Indeed they are proportional to the radion coupling and they must involve M c since in the 
limit of decompactification the radion is part of the higher dimensional graviton, and thus 
is massless. The fact that the radion is typically light and its mass protected was stressed 
in [|l4|] where it was also argued that for M c ~ TeV, the radion would be responsible for 
detectable deviations to Newtonian gravity in present or future experiments. Thus, there is a 



lower bound on M c from short-range gravitational experiments [14, 15, 16]. A radion with a 
mass of 10 -3 eV would induce modification of Newton's law at distances of order 100 /mi (at 
the edge of distances probed by gravitational experiments). No deviation has been measured, 
therefore 

m r > 10~ 3 eV -» M c > 0.8 TeV, (2.12) 

assuming the gluons are bulk fields. The limits for brane gluons are weaker as a result of a 
weaker coupling between brane nucleons and the radion (see below). In these cases, direct 
collider limits will provide a lower bound on M c of order 1 TeV, the precise value depending 
on the extra-dimensional framework [jl7| |. 

2.2 Radion Couplings 



Radion couplings to bulk fields may be simply read from Eq. (2.9), and are characterized by 
the four dimensional Planck scale, mpi. For much of the parameter space, M c < 10 8 GeV, the 
mass of the radion is less than 1 MeV, and decays proceed only into photons (and zero-mode 
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gravitons). These results formally apply only to flat extra dimensions and bulk SM fields. 
Two other possibilities are brane world models 111] in flat extra dimensions, or warped (RS) 



extra dimensions [19|. In the RS case, the radion mass and couplings are both characterized 
by the TeV scale [20. ^lj. This large mass will ruin most of the cosmological applications of 
the radion we are considering, and thus we will not consider the warped case any further (see 
however, Ref |22[|), 

For flat extra dimensions with SM fields on a brane, the radion couplings are substantially 
different because brane fields do not contribute to the fifth components of the stress-energy 
tensor, and thus couplings to radions are induced only through the mixing of 555 with lon- 
gitudinal gravitons [^3|. The net result is that the radion couples only to the trace of the 
brane stress-energy tensor. For gauge fields, this occurs at tree-level proportionally to the 
gauge-field mass, and at one loop proportionally to the (3 function. The brane photon thus 
couples to the radion more weakly than a bulk photon would, with a relative coupling that 
is roughly IQ-k 2 smaller. We parameterize the radion coupling to SM fields by c/mpi, with 
c = 1 for the bulk SM, and c ~ 1CP 3 for the brane SM case. However, zero-mode gravitons, 
as bulk fields, will always couple with full strength to the radion, and thus the brane photon 
case does not substantially affect the radion lifetime, but instead suppresses the branching 
ratio for decay into photons. We will focus on the bulk photon case below, but occasionally 
comment on the brane-world scenario where appropriate. 

2.3 Radion Lifetime 

The decay into two bulk photons implies a radion decay width 

r-r-g (2.13) 

where we have ignored a similar contribution from the decay into zero-mode gravitons. For 
larger M c , decays into other SM fields such as light fermions and hadrons can increase the 
width by a factor of order unity. The radion is effectively stable provided r is larger than 
Hq 1 ~ 4 x 10 17 s, corresponding to 

M c < 7 x 10 8 GeV. (2.14) 

2.4 Radion as the Infiaton 

The radion potential at leading order is V^(0) ~ m 2 (j) 2 /2, for small 4>/mpi. At large cj)/mpi 
this form represents an assumption about the form of the (</>). In order for the radion to 
be the infiaton, the potential energy density of the radion must dominate the energy density 
of the universe. For the radion to be a viable infiaton candidate, it must satisfy the slow-roll 
condition. The quadratic potential is of the chaotic type, so the slow-roll condition does not 
depend on the radion mass, but only on the amplitude of (J). The slow roll condition can be 
stated as a condition on a parameter e: 
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i r\ — 12 i i i i i i i i i i i i i i i i i i 

10 9 10 10 10 11 10 12 10 13 10 14 10 15 

M c (GeV) 

Figure 1: Reheat temperature associated with radion decay as a function of the compactification 
scale. 

which is satisfied for <fi > \[2mpi. As discussed in the next section, this condition turns out 
to be natural for the radion. Note that the requirement 4> > mpi is still consistent with 
a regime where our effective theory is under control. Indeed, the mass of the radion being 
typically much smaller than mpi, the energy density associated with the radion is still much 
smaller than m pl , and therefore we are still away from a regime where 4d quantum gravity 
corrections would dominate. 

One important consideration on inflation is the requirement that the inflaton energy 
density be converted to radiation after inflation and prior to nucleosynthesis and/or the 
freeze-out of dark matter relics. The reheat temperature in a scenario with perturbative 
reheating is given by Trh ~ yJTmpi ~ M 3 / m Pl . We present in Fig. [j] the plot of the decay 
temperature of the radion as a function of M c , (i.e., the reheat temperature of the universe 
in the scenario where the radion is the inflaton). Requiring that the reheat temperature is 
larger than the TeV scale (so that the electroweak candidates for cold dark matter such as 
LSPs and/or LKPs are preserved) leads to the requirement M c > 10 13 GeV. 

The fact that the radion couples gravitationally to any field, and in particular to other 
scalars of the theory, modifies the naive radion scalar potential written above. The effect 
is not negligible if the energy of these scalar fields is significant and leads to an interesting 
inflationary cosmology as we discuss now. 



- 9 - 



3. Radion Mixing and Induced Inflationary Potential 



Let us assume the existence of a scalar field ip other than the radion. If ip is a bulk field, 
according to Eq. (|2.7|) , the total scalar potential is 

7(0, V) = W) + e-*/^"*%ty). (3.1) 

We continue to assume that the form 7,(0) = m 2 4> 2 /2 remains valid at high energies. This 
provides the simplest phenomenological description of stabilized extra dimensions during in- 
flation. Even if the position of the minimum is displaced because of the second term in Eq. 
( |3.1[ ) and thus the size of the extra dimension is modified, V^{4>) ~ m 2 (f) 2 guarantees that the 
VEV of <f> remains finite and is not destabilized during inflation. This potential is plotted as 
a function of <f> for fixed tp and a few choices of r in Fig. |2[ 

When i/j is localized on the 4-d brane, the potential is of a slightly different form: 

V{M) = V^ + Ae-^/^^V^). (3.2) 

In this case the coupling between the radion and ib comes from ,/q T^ u . Here A is an 

1 u ' v induced ^ 

affine function of £, where £ reflects the presence of terms in the (brane) action of the type 

v^wa £ ^ + •••) (3-3) 

which mixes the Ricci scalar with ip 2 . One can suppress the 4> — ip coupling at tree level 
by suitably choosing £, however this coupling will generically reappear at loop level. Conse- 



quently, we expect a potential of the type in Eq. (3.1), and the bulk and brane cases should 
not be qualitatively very different. We will focus on the bulk case below. 

We now investigate the virtues of the potential of Eq. (|3.1|) for inflation. Because of the 



t/j coupling, the VEV of 4> is shifted from zero and satisfies the relation 



^_ P -W/v^ (34) 



\/3m 



where from now on, we refer to eft (and tji) in Planck units. We can easily evaluate ((f)) in the 

case where is a constant; for instance, in the regime where ifi is slowly rolling. We define 
the scale Mj to be V^(ip) = Mf, and in terms of the parameter r defined as 

" = M? (3 ' 5) 

we can express (<p) as 

(4>) = and V(m = SMfr* ^ + ^ . (3.6) 

We plot in Fig. the VEV of cj) as a function of r under these conditions. 
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10" 8 10" 1 10° 10 1 10 8 

0/m pl 

Figure 2: The scalar potential as a function of <f> (in Planck units) for r 
10\ 10°, 1CT 1 , 1CT 2 , 1CT 3 , 1(T 4 (from left to right), r is defined in Eq. ET" 



We can see that in the limit r > 1, the position of the high energy minimum nearly 
coincides with the position at low energy (0 = 0). When r > 1 we have ((/)} ~ l/(3\/3r 4 ) 
leading to V ({(/>)) ps Mf. Therefore, we recover (<p) ~ and the 4> — ip coupling is irrelevant. 
As can be seen in Fig. when r> 1, ^ dominates the potential and plays the role of the 
inflaton. 

In the opposite limit, r < 1, we find ((/>) 3> 1 (see Fig. ||), which leads to a suppression 
of the contribution of ip in the potential. When r <C 1, 4> dominates the potential, V({4>}) ~ 
3Mfr ((j)) 2 /2, and cf> plays the role of the inflaton even though 3> V^. 

A very interesting property of this potential is that even if initially the conditions were 
not satisfied for <f> to be an inflaton, <j> will be driven to be an inflaton through its interactions 
with t/j. If if) carries enough energy density, it inevitably pushes the radion away from (j) = to 
large values of cj>, where the radion can act as an inflaton. Somewhat counter-intuitively, the 
higher the energy density of vp, the more efficiently the radion becomes an inflaton. Because 
of the exponential coupling between <f> and tjj, the higher Mj, the larger the shift of 4>, and 
the larger the eventual suppression of the energy density of tp. Note that during inflation the 
size of the extra dimension is modified. It is larger by a factor e^'' s compared to the low 
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Figure 3: The vacuum expectation value of <j> (in Planck units) at the minimum of its potential as a 
function of r (as defined in Eq. 3.5). The dashed line corresponds to the case where the inflaton lives 
in the bulk. The plain lines correspond to cases where the inflaton is localized in four dimensions with 
suppressed coupling to the radion parametrized by A as defined in Eq. (3^). From top to bottom, the 
solid curves are for A = 0.1, 0.01, 0.001. 



energy case. The mass of the radion is also modified to 

It is informative to compute the value of the ratio m 2 s /H 2 during inflation: 



m 2 eS _ 6 1 + {</>)/ y/3 
H 2 <</>> 2 1 + 2V3/(0) ' 



(3.8) 



In the limit r <C 1, m 2 s /H 2 — > 2y/S/{4>) "C 1. As is well known, in this case the primordial 
fluctuations of 4> are frozen during inflation. This confirms the fact that the radion is a 
viable inflaton candidate when r<l since it can be the origin of cosmological perturbations. 
However, in the opposite limit r > 1, the effective radion mass is too large: m 2 s /H 2 — > 
\/3/(0) » 1. 



- 12 - 



io- 



- I 1 1 1 — I — I — I — I— 



- I 1 1 1 — I — I — I — I— 



-i 1 r~ 




IO" 2 



V e -(0>/V3m pl 



\ 

\ 

\ 



\ 



\ 



\ 



_l I I I I I I L 



J I I I I I I l_ 



_l I l_ 



io- 1 



10° 



Figure 4: Scalar potential (in units of Mf) evaluated at the minimum of <p (solid curve). The dashed 
curve is the contribution from V^> only. The dotted curve is the contribution from ip: V^,e~^l ^ . The 
potential involving the field tp dominates the energy density for r> 1, while for r <C 1, <j> dominates. 



r is defined in Eq. 3.5 



4. Evolution of the Radion Field 
4.1 Evolution During Inflation 

We are now interested in determining the evolution of the radion field in the primordial 
stages of the universe. From the discussion in the previous section, one cannot ignore the 
coupling between eft and if). If we want to follow the cosmological evolution of the radion, it is 
necessary to solve the system of coupled differential equations involving <j> and tp. However, 
for pedagogical purposes, we start with the textbook case and first solve the equation of 
motion for the VEV of (p in the absence of such a coupling : 

<£ + 3F0 + ^(0) = O, (4.1) 

where H = M?/V3m P i is a constant during inflation. In terms of t = Ht and x = m r /H = 
3r 2 , the equation of motion reads 



4> + 30 + x 2 (t) = 0, 



(4.2) 



which has solution 



— 3i/2 

#t) = L_ [(3^0 + 20!) ( e A '/ 2 - e- At "/ 2 ) + A0 O ( e Af/2 + e"^ 2 )! , (4.3) 



2A 



where A = y/9 — 4x 2 , (f>o = 4>(t = 0), <pi = <fr(t = 0). We typically expect <pi ~ x4>o (assuming 
that the initial energy density of <f> is equally distributed between kinetic and potential energy). 
As is well known, for x 3> 1, (/> is exponentially damped during inflation while for x <C 1, <j> 
is frozen and remains unaffected by inflation. However, as we will see now, these conclusions 
are substantially altered when including a coupling between the moduli and the inflaton, even 
as weak gravitational coupling. 

If -0 is a bulk field, the coupled equations to be solved are 



k'2 



2 ' M^j + ^ + ^"» 



pi 



^ + v^j Y + T +c rT + -^-j +d ^+ — — 2 =0, 

where eft and ip are in Planck units, we have defined a dimensionless time variable t in terms 
of a convenient mass scale M as t = M.t, and we have introduced the constants 

Cr = —tk , d r = -f- , di = -fj. (4.5) 
M 2 MM 

Note that the addition of the decay rate terms in the equations of motion is a an effective 



description which is valid only when oscillations of the fields have started [ 24 1 . 

We define D to be the value of (j) in Planck units when the field oscillates and its energy 
density starts to decrease as a nonrelativistic component: 

D = MhA. (4.6) 
mpi 

The factor D may be interpreted as the damping factor if the initial value of the radion is 
Planckian. Our goal is to evaluate by which amount the radion field can be damped during 
inflation for a given value of r, so that we can investigate the resulting radion energy density 
today. 

The damping of the amplitude of generally depends on the form of the potential for ijj. 
However, some general statements can be made. Indeed, what typically happens can be easily 
understood by inspection of Fig. ||: Because of the coupling between (j> and tp, the position 
of the minimum for <p is displaced from its low energy position. During inflation, eft is pushed 
to its high-energy minimum and we naturally expect that once inflation ends, oscillations in 
(f) start with an amplitude of the order of the VEV of <fi as given in Fig. ||. We checked that 
this is indeed what typically happens using different types of potentials for tp. 
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r=0.5 




r=10 




Figure 5: Examples of the time evolution of 4> and ip using the potential for ip of the form V^(ip) = 
m^ip /2. In these specific examples, r = 0.5 in the top panel and r = 10 in the bottom panel 
(r = \/ m r I rrhp ) . Initial conditions were chosen such that ^0 = 5, and tfio to satisfy the equipartition 
condition. 



As an illustration, we show in Fig. || the time evolution for the <fi and ip fields in the case 
where the potential for if) is of the large-field (chaotic inflation) type. 2 What is shown in the 

2 An interesting possibility for the origin of the "precursor" inflaton ip is A$, the higher component of a 
gauge field §§]. 
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Figure 6: D(r) as derived in the chaotic case. For initial conditions we have assumed i/jq = 5, and 
have chosen <^>o to satisfy the equipartition condition. 

figures are the values of (c/> 2 ) 1//2 and (ip 2 ) 1 / 2 : 

{4> 2 ) 1 ' 2 = y/p*/™? (4.7) 
(V 2 ) 1/2 = ^p^/mlexp(-<p/V3), (4.8) 

where the energy densities capture both the kinetic and potential energies. During inflation, 
((J) 2 ) 1 / 2 = <f>/y/2 and (ip 2 ) 1 / 2 = Ve"* /2v ^/2. After inflation, when the fields oscillate, the 
averages defined above correspond to cycle averages of the field. 

Note that in the case of the chaotic potential rn^ip 2 /2, r ~ ^Jm^jm^,. For r to be larger 
than unity would require that the mass of ip be somehow less than the radion itself. This 
may be difficult to achieve naturally, however, we consider this possibility since we find that 
the chaotic potential for ip provides the largest damping during inflation. (See the D(r) plot 
derived in this case in Fig. ||). As initial conditions, one may choose equipartition of the 
potential energy densities between <j) and vp, leading to 

e-^ = rUl/^l (4.9) 
and the initial condition for ipQ has to be chosen such that tpo satisfies a slow-roll condition. 
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In the following sections we will express cosmological constraints as general bounds on 
D, and will illustrate how these constraints typically translate into r, remembering that for 
r < 1, we expect D will naturally be of order unity. 

4.2 Validity of 4D effective field theory 

Throughout the discussion we have assumed an effective 4-dimensional description. The 4d 
effective field theory is an approximation which neglects the effects of higher KK modes. 
Those can be excited at high temperature. For the 4d description to be valid, we should have 
the temperature smaller than M c . This requirement is in addition to the 5d cut-off A which 
characterizes the scale at which the 5d gauge theory (including KK modes) becomes strongly 
coupled, and loses predictivity. 

The definition of the frontier between the 4d and 5d descriptions in the inflationary phase, 
where the universe is dominated by vacuum energy, is somewhat subtle. During inflation it 
is possible to define a temperature. In the exponentially expanding universe particles are 
created with typical momentum of the order of H. In this phase, the criterion for the validity 
of the 4d effective description is H < , where Lp^j/s = e^/^* /M c . This leads to a 

constraint on M c , shown in Fig. [?]. The corresponding constraint is weak, and is only relevant 
for small values of r, which, as we will see in the next section, will be disfavored anyway 
because of the problem of the overclosure of the universe by radions. 

The reheat temperature satisfies Trh < Mj where the bound is saturated in the case of 
preheating. We can define a parameter V = Trh /Mj which measures the efficiency of the 
transfer between the inflaton energy to radiation energy. Requiring an effective 4d description 
(Trh ^ M c ) leads to V < r. Note that in the case of perturbative reheating Trh ~ X^/mjmpi 
where A is the inflaton coupling to normal matter and mj is the inflaton mass. The require- 
ment that the theory appear four dimensional thus translates into mj < r 2 A~ 2 10 11 GeV. 

5. Radion as Dark Matter 
5.1 Coherent Production 

We now turn to the evolution of the radion field once inflation is completed. The story 
begins with the energy density of the radion after inflation (determined by the parameter D 
discussed in the previous section). From the time oscillations commence until today, 3 the 
radion energy density p r redshifts as nonrelativistic matter, so that 

a 3 

Pr(t ) = Pr(tosc)^f, (5.1) 
«0 

where the subscript '0' refers to the present, and 

Pr(tosc) = \m 2 r 4> 2 {tosc) = \d 2 M*. (5.2) 

3 If the radion is dark matter it must have a lifetime longer than the age of the universe, which requires 
M c < 7 x 10 8 GeV. 
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r 

Figure 7: For the Hubble radius to be larger than the size of the extra dimension Lphys — 

during inflation, M c must be in the lower region limited by the dashed line. This does not lead to 

strong constraints in the preferred region r > 1. 

After inflation, until reheating (or preheating), the energy densities in both the inflaton 
and the radion decrease as a nonrelativistic component. We can express p r (to) as 

Pr (t ) = |d 2 M c 4 %%. (5-3) 
i a RH a 

Before reheating, a 3 oc t 2 . Since t 2 ~ H~ 2 oc p 1 , we can replace the factor o^ sc /a RH by 
Prh/Posc- The value of p RH is simply (Tr 2 g^(T RH )/30)T^ H . The value of p osc is not so 
straightforward. 

Here we will use the results for the chaotic inflation example of the previous section. If 
r < 1, then the radion is the inflaton. In this case inflation ends when (f> ~ mpi and p osc ~ 
m 2 m 2 Pl . Also in this case, D ~ 10 -1 . If r > 1 the energy during inflation is approximately 
Mf and D ~ 10 -2 r~ 4 (as can be seen from Fig. |6|). Using these values, along with the fact 
that a RH /al = g*(T )T$ / 'g^(T RH )T^ H , we find 

^(to) = ^(To)T 3 r^|^ r °; 4 3 r r ~[ (5.4) 
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Using T = 2.4 x 1(T 13 GeV, g(T ) = 3.36, and p c = 3.9 x 1(T 47 GeV 4 , this leads to a value 
of £l r of 

' 2 x 10 6 (T«ff/GeV) r < 1 

6 x 10 4 (Tjuj/GeV) r -4 r > 1. 



If the radion is stable, we require f2 r < 0.3, and the inequality is saturated if the radion 
is dark matter. The limit to Q r results in the limit Trh < 5 x 10 _7 GeV for r < 1, and 
Trh < 2 x 10 _5 r 4 GeV for r > 1. Clearly the r < 1 case is ruled out. The parameter r must 
be sufficiently large to avoid overclosure: 

r>15(|f) 1/4 . (5.6) 

We now discuss some potential observational effects associated with the oscillation of 
the radion 4 . Coherent oscillations of the radion induce variations of gauge couplings (in a 



way described in Ref. [14], for instance). For this effect not be observable, we require the 
oscillation period m~ l to be much smaller than the Hubble time. There are constraints 
on the variation of the gauge couplings at the time of nucleosynthesis. Applying the above 
requirement to the nucleosynthesis epoch leads to H nuc [ eo <C m r , i.e., T nuc i eo <C M c , which 
is always satisfied. Note that as soon as structures start to form, coherent oscillations will 
get distracted by the development of gravitational potentials. Therefore, radion oscillations 
would no longer be coherent today and would not be seen as variations of the couplings, but 
instead, as nonrelativistic matter. 

5.2 Thermal Production 

Now we turn to the thermal production of radions at reheating. The number density of 
radions produced thermally is obtained by solving the Boltzmann equation 

(in 

-± + 3Hn r = -{av)n 7 (n r - n e r q ) - T(n r - n e r q ), (5.7) 

where re 7 represents the number density of a typical light species in the reaction radion + 7 — > 
everything. The cross section is expected to be approximately m™, and the decay width is of 
order m^/m^. 5 The feeble radion interactions with matter imply that the radion has great 
difficulty reaching thermal equilibrium. We have also seen above that the coherent production 
of radions need be very small. Thus, initially, n <C n eq . In terms of Y = n r /s where s is the 
entropy density, and x = m r /T, the Boltzmann equation becomes 

dY m r (C A \ 

— = —5- + xC D Y eq , (5.8) 

ax mpi \x z J 

where Ca and Co are dimensionless numbers. 



4 See also jjj, g 

5 There are potentially large dimensionless numbers in the cross section and decay width; e.g., the factor of 



1927T appearing in the decay width given in Eq. (2.13) 
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If the radions are relativistic (x <C 1), Y eq ~ x , and 

Y ~ — . (5.9) 
mpi 

We note again that there is a potentially small (say 1CP 2 ) dimensionless number on the 
right-hand-side of the above equation. 

We now examine the case where the radions are nonrelativistic at the time of reheating 
(x S> 1) so that Y e q ~ x 3//2 exp(— x) and we can ignore the first term in the right hand side of 
Eq. (|Q|). In this case 

in 

Y — I(xrh), (5.10) 

mpi 



with 



oo 



T{x RH )= / x 5 / 2 e~ x dx, (5.11) 

where xrh = m r/TpH is x at the time of reheating. 

We find that f2 r from thermal radions (for relativistic radion production) is 



-Trh ( M c 



2 



f] r = 10 9 -^ —1 , (5.12) 
GeV \mpi/ 



while the value for nonrelativistic radion production is 



The limit Vt r < 0.3 becomes 

2 



M 4 

n r = lO^eV- 1 -^ x I{ XRH ). (5.13) 



/ M c YTrh 
UxlO^Gev J GW <1 (5 " 14) 

for relativistic radion production, and for nonrelativistic radion production 

4xl0n G ev ) 4 ^' < L < 5 ' 15) 

This last bound is always satisfied, while the bound from relativistic radions only leads to a 
very weak constraint on the reheat temperature. Therefore, there is no consequence, as far 
as dark matter is concerned, from thermally produced radions. 

Radionactive dark matter would be challenging to detect directly, though it may have 
implications either for gamma-ray spectra (through the eventual decay of the radions) or 
isocurvature perturbations. 
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6. Radions and Cosmological Perturbations 



6.1 Radion Isocurvature Perturbations 

A possible effect from the stable radion is on the CMB, an issue we now investigate. Indeed, 
in a way reminiscent of the axion, there are isocurvature perturbations associated with the 
radion as dark matter. To discuss this, we begin with a few definitions. The curvature 
perturbation associated with an energy density component pi is defined by 

Ci = -H 6 -^. (6.1) 

Pi 

For matter, £ m = \5p m /p m and for radiation £ 7 = ^Sp^/p^. The adiabatic mode of the 
perturbation is defined by Q = £~. The nonadiabatic or isocurvature modes arise from local 
inhomogeneities in the equation of state of the energy density. They may be specified relative 
to £ 7 as 

*^-^ = 3fc-Cr) (6.2) 
Pi 4 p 7 

The large-scale temperature anisotropy can be related to the primordial adiabatic and isocur- 
vature perturbations by (see for instance section II of Ref. |27|): 

5T . . 1 ^ad 2 ( Pc c, , Pb r\ , 1 Pv 



W - -^S c + ^S b \+ -^-^—Sr (6.3) 

J 5 5 \p m p m J 15 p 7 + Pu 

where (, ad measures the primordial adiabatic perturbation, S c , S b and S u the isocurvature 
perturbations in cold dark matter, baryons and neutrinos respectively. Here p m is the total 
matter energy density, p m = Pc + Pb- 

We are interested in constraining S c . From an observational point of view, what is 
constrained is the quantity: 

ceff 

B = P, (6-4) 

where 

Sf = S b + ^S c (6.5) 
Pb 



is the effective nonrelativistic term. We use the constraints derived at 95% C.L. in Ref. [27] 

-0.46 < B < 0.35. (6.6) 
In the case where the radion is the dark matter we have 

"Sc = S ra dion = ^{Cradion Cy) = = & j (6-7) 

P<t> 

where, at first order, 

§ iso = 
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If the radion is practically free during inflation, with small effective mass m 2 <C H 2 , it follows 
that on superhorizon scales 

6<t) H* . . 

(6.9) 



(j> 2lT(j)* 

where and are H and <f> evaluated at the time of horizon exit, and 

giso = I±= ^* ^ ( g 1Q) 

tt4>* \/3'KD*m i p l 

where D* is 4>*/mpi. Equation fl6.6|) leads to the restriction \5 tso \ < 10 _2 C 7 , which leads to 

y/V* < 10~ 5 D*m 2 Pl . (6.11) 

This constraint, while interesting in its own right, is not relevant for our current picture. As 
we will see in the next section, m 2 /H 2 < 1 during inflation only if r < 0.4, which typically 
results in D of order unity. But if D is of order unity, the limit in Eq. ( |6.11| ) is weaker than 
the limit on V* from the limit to the gravitational wave background produced in inflation. 
However, we note that there may exist unconventional inflationary mechanisms where this 
limit is useful. 

To summarize, we have seen that when the compactification scale is below 7 x 10 8 GeV 
and the radion is effectively stable, to avoid overclosure of the universe by radions we must 
have r > 15(Trh /GeV) 1 / 4 . 

If D and r are properly correlated, radions could in fact play the role of cold dark matter, 
though in the models of inflation we have considered this would require fine tuning. 

6.2 Radion as the Curvaton 

It is now widely believed that the origin of structure in the universe is a primordial per- 
turbation already in existence when cosmological scales start to enter the horizon. In the 
standard picture for the generation of cosmological perturbations, the curvature perturbation 
C, is generated during inflation through the perturbation of a single component inflaton. 



In the curvaton scenario [11] the curvature perturbation produced during inflation is 
initially negligible, but the curvaton field acquires "isocurvature" perturbations during in- 
flation, which is converted to curvature perturbations when the curvaton density becomes a 
significant fraction of the total density. In this case, £ is given by 

C = fo x — (6.12) 

where fo is the fraction of the total energy density (radion r and radiation R) of the universe 
in the radion-curvaton when it decays: 



Pr 



Pr + %PR 



(6.13) 
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When the radion decays, it reheats the universe and imprints its perturbation into photons 
and matter. The degree of nonadiabaticity is given by the fraction fjj which reflects how much 
of the total energy density is carried by the radion at the time of its decay. If this fraction 
fo is unity, there will be no residual isocurvature perturbations. On the other hand, if the 
radion does not completely dominate the energy density of the universe at the time of decay, 
there will be residual isocurvature perturbations. 

If m 2 <C H 2 during inflation then the spectrum of the fractional field perturbation is 

K 5: <fU4) 

and that of the fractional energy density perturbation is 

VlH^E^^S™ . (6.15) 

TT(ft* 

It follows that the prediction of the radion-curvaton model for the spectrum of the curvature 
perturbation is 

1/2 ^ fD 6lso ^ [dH^ (g 16) 

<> 3 

This is to be compared with the COBE measurement of the CMB quadrupole anisotropy 

V 1 ^ 2 = 4.8 x 1(T 5 . (6.17) 

It is now clear that the radion cannot easily play the role of a curvaton. Indeed, from 
Fig. [8|, we see that the effective mass of the radion is smaller than H only for r < 0.4. On 
the other hand, Fig. ^| indicates that for small r, (ft dominates the total energy density during 
inflation, in conflict with the curvaton scenario. In fact, this leads back to the model of the 



radion as the inflaton (and providing the primordial fluctuations) discussed in Sec. 2.4. 

It is possible that ift could play the role of the curvaton. This would require that tft decays 
after (ft, which may be unlikely for low compactifications scales since generically a scalar field 
would have gravitational interactions, and thus at most a lifetime of the same order as the 
radion. This would require Xm^p < M^/rrip,, where A characterizes the tft interaction strength. 

Let us now see if the radion can play the role of a modulon. 

6.3 Radion as the Modulon 

As mentioned in the introduction, another alternative to the standard picture of generation 



of cosmological fluctuations was recently proposed [12], which also makes use of a light scalar 
field that controls the coupling of the inflaton to matter. To distinguish it from the curvaton 
scenario, we refer to the field responsible for the generation of fluctuations in this mechanism 
as a modulon. In the modulon scenario, the coupling A between the inflaton and ordinary 
matter is not constant but is given by the fluctuations around the VEV of the modulon. If the 
modulon is light enough, it acquires primordial fluctuations during inflation of order H*/2tt. 
These spatial fluctuations will then correspond to spatial fluctations in the coupling A. In the 
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Figure 8: The ratio m^ s /H 2 during inflation (evaluated at the high-energy minimum of the potential 
for cf>) . Fluctuations in (f> are frozen and can be at the origin of cosmological perturbations only for 
r < 0.4. 



following we will denote the modulon by 4> since we want to check whether the radion can be 
a modulon candidate. The (local) reheat temperature depends on this coupling according to 

where Fj ~ X 2 mj is the decay rate of the inflaton and mj its mass. In the modulon scenario, 
where A is given by 

A(</>) = A (l + + ...). (6.20) 

As the radion controls the size of the extra dimension, it controls the inflaton coupling to 
matter as in Eq. fl6.20D , where Ao is related to a dimensionful, 5d parameter by Ao = A5/L. 
Small variations in the radion VEV thus indeed lead to variations in the coupling between 
matter and the inflaton, with M = y/Smpi. The fraction of the coupling controlled by the 
fluctuating radion is (j)*/mpi where 0* = D^mpi is the value of <fi at the time of horizon exit 
so that Eq. ( |6.19 ) can be written as 

5T 0* H* H* 



T mpi27T(j)^ 2Trmpi 



(6.21) 
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a contribution typically much too small to account for observed inhomogeneities. Thus, in 
models with flat extra dimensions, the radion cannot effectively play the role of a modulon. 
For this mechanism to be effective, the radion coupling should be much larger, as for example 
in warped compactifications. In the Randall- Sundrum geometry, we would arrive at ST /T ~ 
H*/l TeV. In order to explain the COBE result, this would require an inflation scale around 
10 s GeV, which is neither the fundamental UV nor the IR scale of the Randall-Sundrum 
model, but could be generated by an intermediate brane . 



7. Late Decay of Radions and KK Gravitons 

The decay of massive relics after Big Bang Nucleosynthesis produces energetic photons which 
may dissociate light elements and modify the predictions of their abundances. The constraint 
on the decay of heavy relics from BBN is (see for instance ) 

mY < 1(T 12 GeV (7.1) 

where Y = nx / s is the ratio we would have observed today if particle had not decayed, and 
s ~ T 3 is the entropy per comoving volume. This bound literally applies for particles whose 
lifetime is 10 8 s; bounds for larger lifetimes are somewhat relaxed. Thus, by taking 10~ 12 



GeV on the right hand side of Eq. (71) we are making a conservative estimate for the allowed 
parameter space. Note that the constraint Eq. ( |7.l[) applies only if the decay takes place after 
nucleosynthesis . 

7.1 Late Decay of Radions 



From our estimate of the radion lifetime, Eq. ( 2.13 ), we find that the radion decays after 
nucleosynthesis (r > 10~ 2 ; T < 10 -22 GeV) provided 

7 x 10 8 GeV < M c < 10 12 GeV, (7.2) 

where the lower limit on M c comes from the requirement that the radion lifetime is smaller 
than the age of the universe. In order to explore the implications on radion cosmology , we 
now estimate Y r = n r /s at the time of reheating (which is the same as the ratio we would 
have observed today if the expansion is adiabatic) . 
The bound in Eq. ( |7.ip may be written as 

nj. = P L = Qr Pc^ 10 -i2 QeV _ (73) 
s s s 

Using p c /s ~ 10 -9 , we obtain the constraint fl r < 10~ 4 . 

For the chaotic inflationary potential, coherent radion production [using Eq. ( |5,5D 1 yields 



T RH < f 5 x HT 11 r < 1 
GeV ~ I 2 x lCr 9 r 4 r > 1. 



(7.4) 
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Recall that this bound only applies in the range 7 x 10 8 GeV < M c < 10 12 GeV, so that the 
radion decays after nucleosynthesis. Obviously the case r < 1 is again excluded. The bound 
on Trh for r > 1 is plotted on figure ||. 

The constraint in Eq. fl7,l| ) also applies to radions of thermal origin: 

from the thermal production of relativistic radions in the radiation era. If M c ~ 10 9 GeV, 
then Trh < 10 6 GeV. If M c ~ 10 12 GeV, then the bound becomes T RH < 1 GeV. 

For nonrelativistic radion production [using Eq. (pTIaD l the bound in Eq. ([O]) becomes 

^„» £ (i^) 4 

This bound imposes constraints on T R h for the larger region of M c . Below M c < 3 x 10 10 
GeV, there is no bound on T RH . At M c ~ 10 12 GeV, the bound is T RH < 35 TeV, for 
M c ~ 10 11 GeV, the bound is T RH < 800 GeV and for M c ~ 5 x 10 10 GeV, it is T RH < 320 
GeV. 

7.2 Late Decay of KK gravitons 

KK gravitons (G 1 ) also couple very weakly, ~ rripj. They have a mass characterized by the 
size of the extra dimension, vtlkk ~ M c and therefore a decay width of the order T ~ /m Pl , 
though the specific details of the expression depend on whether KK parity is well conserved, 
which would forbid decays directly into zero modes of SM fields. Continuing with dimensional 
analysis, the decay is after nucleosynthesis if M c <8x 10 4 GeV. On the other hand, we only 
consider values of M c above several hundred GeV, as lower values are generally excluded by 



collider constraints [17]. Thus, constraints on KK gravitons from BBN will only apply in the 
narrow range 

300 GeV < M c < 8 x 10 4 GeV (7.7) 

In this range of values for M c , the radion is stable, and to avoid overclosure of the universe 
by radions, r (D) and T R h must be chosen appropriately in Eq. ( |5.5| ). 

We assume the universe is radiation dominated at reheating, and compute bounds from 
the thermal production of KK gravitons. As for the radion, (av) ~ and many of the 

radion results can be carried over with m r — > M c . We begin with the case of relativistic 
KK gravitons, T RR > M c . In this case, thermal production of KK gravitons leads to Y ~ 
Tftn/mpi and the BBN bound results in, 



^10 6 GeV\ 



Trh < ( , r ) GeV (7.8) 



from the late decay of thermally produced, relativistic KK gravitons. For M c at the TeV scale, 
this roughly constrains T RR to be of order M c . For somewhat smaller extra dimensions, this 
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Figure 9: The constraints from overclosure by radions (produced non-thermally) -dashed curve- as 
well as from BBN -solid curve- due to the late decay of radions, translated into a bound on Trh as a 
function of r. The bound from overclosure (see Eq. 5.5) applies in the range M c ~ 10 3 — 7 x 10 s GeV 
and the one from BBN (Eq. O) in the range M c ~ 7 x 10 8 - 10 12 GeV. 



requires a smaller reheat temperature, and indicates that KK gravitons cannot in fact be 
relativistic for M c > 1 TeV. 

For nonrelativistic KK gravitons, Y ~ M c /mpiXRH^, where now xrh = M c /Trh. Thus, 
the bound on the reheat temperature is, 

which for all M c of interest results in a bound of Tru < M c /afew 6 , consistent with the 
possibility of KK dark matter ||, and in fact with the general WIMP hypothesis, for which 



the freeze-out temperature is typically a few tens of GeV. Note also that Ref. [31] studied the 
case where KK dark matter was made of KK gravitons rather than KK hypercharge gauge 
bosons, with the relic abundance determined initially by freeze-out of B . They derived 
bounds from BBN on this scenario coming from the eventual late decay of B 1 into G 1 and a 
photon. 



6 This point was realized some time ago by Abel and Sarkar 
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To finish this section, we discuss the possibility of nonthermal production of KK gravitons. 



As stressed in [32|, the gravitational production of particles during inflation can be very 
efficient and can lead to much stronger bounds on the reheat temperature than the ones from 
thermal production. However, for this mechanism to be operational, we need the mass of the 
particle to be smaller than H during inflation. In our case, this corresponds to e~^' 3 M c < H. 

From Fig ^, this leads to a bound on M c only for very small values of r. On the other 
hand, because of the bound from overclosure of the universe we typically have r > 1 so that 
no KK gravitons are expected to be produced nonthermally. Note also that from Fig. ||, 
there is essentially no gravitational production of radions during inflation for r > 1. In any 
case, the bound from overclosure of the universe by radions is not the only reason why a 
small hierarchy between M c and Mj is required. A large enough reheat temperature, to have 
WIMP dark matter, also requires large values of r. 

7.3 Diffuse Gamma-ray Signal 

Finally, for a very long lived radion, there might be some interesting observable diffuse gamma- 



ray signal 2£, 33]. This applies for life times larger than 10 14 s, i.e., for 

M c < 3 x 10 9 GeV. (7.10) 

On the other hand, the radion becomes stable for r > 10 17 s, i.e., for M c < 9 x 10 8 GeV. The 
observable signal therefore corresponds to a tiny window of values for M c . This is because the 
decay rate depends on the sixth power of M c . Note that there is no similar expected signal 
from decaying KK gravitons (however see |3l| for the special case of UED) since r > 10 17 s 
corresponds to the disallowed region M c < 0.4 GeV. 

8. Conclusions 

In this work we have studied the cosmological properties of the radion field as it appears in 
models with flat or almost flat extra dimensions. We find that the potential roles the radion 
may play in cosmology are very sensitive to the scale (and model) of inflation. Thus, our 
work explores the challenges involved in building a realistic cosmological model for theories 



with flat extra dimensions. Figure 10 summarizes our results. 

Its origin as a component of the higher dimensional graviton tensor dictates the radion's 
couplings. Similarly, the mass of the radion, in the absence of fine tuning, is generally related 
to the compactification scale M c by m r ~ Ml/mpi. Thus, once one has chosen the size of the 
extra dimension, the natural radion properties are determined. Of course, it is possible that 
the radion potential is itself fine tuned by some mechanism we currently do not understand. 
In that case, its couplings to other fields may still be of gravitational strength, but its mass 
could possibly differ from our estimates. One avenue of future research would be to explore the 
compactification dynamics that could naturally make the radion, for example, more massive 
than the naive estimates, such that it decays at early times and evades some cosmological 
constraints. 
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Figure 10: A summary of the cosmological significance in various regions of radion model parameter 
space. M c is the compactification scale and Mj is the scale of the potential of the scalar field responsible 
for inflation. We note that the application of effective field theory breaks down in the region Mi > 
10 3 Mc, as well as the region where the Hubble radius is smaller than the size of the extra dimension. 



We have written the scale of the radion potential as M c and used a naturalness argument 
to identify L ~ M^T 1 , without specifying the compactification dynamics. It is interesting 
to compare this with existing stabilization models, and the radion masses which result. In 
Ref. [34], it was shown that casimir effects will stabilize the extra dimension (and in fact, 
are always present for a compact space). This leads to a potential which is natural in our 



sense, and an estimate of the radion mass similar to Eq. ( 2.11 ) |14| . An alternative (in 
6 dimensions) is to invoke a bulk cosmological constant and a set of gravitating 3-branes, 
balanced against quantum effects from KK mode gravitons |35|] . This also leads naturally to 
the relation between the scale of the potential (characterized, i.e., by the bulk cosmological 
constant) and L which we have assumed. Unlike the TeV -1 (or smaller) extra dimensions 
we consider, the aim of those works was to generate large (mm) size extra dimensions, for 
which the radion is too light, in conflict with precision measurements of gravity. In order to 
generate acceptable large extra-dimension phenomenology, one may introduce a bulk £7(1) 
gauge theory and stabilize with its trapped magnetic flux [35, 36 1. An additional parameter 
allows one to tune the radion potential such that the simple relation between the radion 
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mass and M c is disturbed, with the free parameter given by the U(l) coupling. Another 
way to avoid the natural relation m r ~ M^/mpi was suggested in |37[| where the stabilizing 
potential depends only logarithmically on the radion field so that a large hierarchy between 
M c and R^ 1 can be generated without fine-tuning. The idea works for an even number 
of dimensions only since it relies on sets of two transverse dimensions (for the logarithmic 
dependence). In addition, bulk supersymmetry is needed to protect against the generation 
of a bulk cosmological constant which would induce power law corrections to the effective 
potential for the radion and spoil the picture with logarithmic potentials. 

Concerning the radion coupling to an inflaton (or any scalar field), the result is an in- 
teresting inflationary dynamics through which the radion may be pushed to dominate the 
energy density of the universe and play the role of the inflaton. Given the scale Mj which 
characterizes the potential of another scalar field, the cosmological evolution of the universe 
depends crucially on the ratio r = M c /Mi where M c is the compactification scale. The cou- 
pling to the inflaton will generally push the radion VEV away from its low energy minimum, 
in way strongly reminiscent of the coherent production mechanism of the axion. Models with 
low compactification scales (less than about 10 9 GeV) possess radions which are stable on 
cosmological time scales, for which there are strong overclosure constraints. When inflation 
ends, energy density is stored in radion oscillations which can easily dominate the energy 
density of the universe at late times. This strong moduli problem requires that, for models 
of conventional inflation, the inflation scale must be below the compactification scale r > 1. 

An alternative is the possibility that the radion could account for the dark matter required 
by cosmological measurements. Radion dark matter is probably challenging to detect directly, 
though evidence might be available from the gamma-ray sprectrum (through the eventual, 
rare radion decay) or from isocurvature perturbations in the CMB. 

We have examined the possibility that the radion could play the role of either the curvaton 
or modulon, and find that generally it is not suitable. It typically does not have large enough 
density fluctuations to be either a curvaton (unless it is already the inflaton) or a modulon. 
Finally, we have explored the possibility that the radion could decay after nucleosynthesis, 
and examined the stringent bounds which result. These bounds require both a much smaller 
abundance from coherent production and require a low reheat temperature to avoid thermal 
production mechanisms. 

8.1 TeV Extra Dimensions 

Models with TeV -1 extra dimensions have attracted much attention lately. Model builders 
in particle theory have made use of TeV -1 extra dimensions to suggest new mechanisms 
for breaking symmetries, naturally generating small numbers, and addressing the hierarchy 
problem [§, ||, ||| using orbifold boundary conditions. The natural scale for which these 
mechanisms to operate is the natural scale of physics beyond the SM. Thus, compactification 
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should be at the TeV scale 7 . Thus far, little has been said 8 about the cosmology of these 
models, aside from the possibility of Kaluza-Klein dark matter, which is only an option in 
the special case of universal extra dimensions. 

In the present work, we have derived additional and stronger constraints from the radion 
physics. We showed that TeV -1 extra dimensions generally have an overclosure problem, 
because the radion is effectively stable and easily dominates the energy density of the universe 
at late times. The radion energy density can be sufficiently damped provided the inflation 
scale Mj is at a TeV. The radion being stable, reheating would have to come from a different 
field which would imprint its fluctuations into the radiation. 

To summarize, in models with TeV -1 extra dimensions: 

• The radion is effectively stable. 

• To avoid overclosure of the universe by radions, either the inflationary model must be 
rethought, or the inflation scale must be below the TeV scale i.e., r > 1. In this case, 

is subdominant during slow-roll and (j> cannot be the inflaton. There is a narrow 
window of parameter space in which the radion can account for the observed abundance 
of dark matter. 

• Inflation model building is challenging: A reheat temperature of order 100 GeV is 
preferred for baryogenesis and also for dark matter (unless the radion is dark matter — 
the LSP or LKP dark matter candidates actually require Trh ~ 100 GeV). For Trh ~ 
TeV~ Mj, this necessitates a very efficient reheating, and almost certainly preheating. 

In the last years, particle physics model building has moved into theories with TeV extra 
dimensions. We have explored the cosmological challenges associated with these models, and 
by seeing where constraints arise, define the properties that a solution of these challenges 
should have. One way out of these constraints points towards non-toroidal compactifications, 



like, for instance, warped geometries [19, 20, 21] or internal compact hyperbolic spaces for 



which m r ~ R 1 |39[| . The future of extra-dimensional cosmology is both challenging and 
exciting. 
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